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Abstract 

A study of the effects of non-extensivity on the modehing of atomic 
physics in hot dense plasmas is proposed within Tsallis' statistics. The 
electronic structure of the plasma is calculated through an average- 
atom model based on the minimization of the non-extensive free energy. 

1 Introduction 

For astrophysical applications, as well as for modelling laser-produced plas- 
mas, equation-of-state data are required over a wide range of physical con- 
ditions. When the ions are strongly coupled {i.e. when Coulomb interac- 
tion between ions becomes larger than the thermal kinetic energy), and the 
electrons degenerate {ksT < ep, where ep is the Fermi energy, T the tem- 
perature and ks the Boltzmann constant), the plasmas combine features 
of hot matter, such as ionization and fluid behaviour, and characteristics 
of cold matter, such as electron degeneracy. Some electrons are bound to 
the nucleus, and others are free (or delocalized) . The delocalized states are 
also sometimes called "scattering states" . The electrons can be treated with 
various degrees of complexity. The idea of the average-atom model is to 
consider only the mean configuration of the atom. In this picture, the atom 
is viewed as confined to the Wigner-Seitz (WS) sphere, which is immersed in 
a homogeneous jellium of delocalized electrons neutralized by a continuous 
background. In the approach proposed by Rozsnyai [1], the bound elec- 
trons are treated quantum-mechanically, and the delocalized electrons as a 
Thomas- Fermi fluid p] . This simple but efficient model has been used in nu- 
merous approaches of plasma structural and radiative properties. The first 
purely quantum approach to an average-atom model is due to Liberman [3] . 
In this model, the WS sphere appears explicitely as a cavity into which non- 
central ions can not enter. The delocalized electrons are treated quantum- 
mechanically, which results in Friedel-type oscillations of the self-consistent 
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electron density and a potential extending beyond the WS radius. Recently, 
a new version of that approach, including efficient numerical methods, was 
developed by Wilson et al. [1]. A great progress in the field was brought by 
R. Piron and T. Blenski [5], who developed a numerical code which provides 
the average-atom structure and the mean ionization self-consistently from 
variational equations. This work enables one to clarify the thermodynamic 
consistency issues in the existing average-atom models. For a sake of sim- 
plicity, in the present work we consider an average-atom model similar to the 
one proposed by Rozsnyai, in which bound electrons are treated quantum- 
mechanically (by Schrodinger equation with relativistic Pauli corrections 
[6]), and delocalized electrons semi-classically (Thomas- Fermi model) |7]. 

In standard thermodynamics of systems in local thermodynamic equi- 
librium (LTE), quantities such as energy and entropy are extensive, which 
means that they are proportional to the size of the system. The lack of 
adequacy of Boltzmann-Gibbs (BG) entropy is related to the breakdown 
of the extensivity. More precisely, BG statistics fails for instance when 
a system includes long-range interactions [51 [9], long-time memory effects 
(non-Markovian process) and/or evolves in a (multi)-fractal space. Tsallis' 
form has been widely used in various fields of physics and is considered as a 
possible framework to deal with non-extensive settings. Lima et al. \\{)\ [TT] 
have studied the kinetic foundations of Tsallis' statistics through a modified 
Boltzmann transport equation satisfying a modified i/theorem. 

When a physical system is made of a large number of identical sub- 
systems, the thermodynamic identity 

{E) -TS- fi{N) = -PV (1) 

is verified, quantities (E) (average energy), (A^) (average number of par- 
ticles) and S (entropy) are extensive and the grand potential provides a 
simple way to calculate the thermodynamic quantities. On the other hand, 
in the case of average-atom models, equation ([1]) is a priori not necessarily 
verifiecil since only one nucleus is taken into account and (N) must be equal 
to the atomic number Z. In the framework of the average-atom model, 
the system can not be understood as a statistical sum of many identical 
subsystems. 

In the present work, we propose to investigate how non-extensivity af- 
fects atomic-structure calculations in hot and dense plasmas within Tsallis' 
statistics. To our knowledge, it is the first time that non-extensivity and 

^It is also the case for the Virial theorem: 

3PV = 2Ek + Ep (2) 

which requires extensive quantities, and therefore a system made of a large number of 
identical sub-systems. In Eq. ((2]), P is the pressure, V the volume, Ek the kinetic energy 
and Ep the potential energy. 
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Figure 1: Comparison between Fermi-Dirac distribution and the non- 
extensive orbital occupancies fq for different values of parameter q. 



quantum shell-structure effects are taken into account simultaneously in a 
self-consistent procedure. In section [21 Tsallis' entropy is introduced [12j . 
In section [3l the self-consistent calculation of the electronic structure (ener- 
gies, occupancies and wavefunctions of the orbitals, etc.) is presented. The 
impact of non-extensivity on ionic distributions is illustrated in section HI 
Section [5] is the conclusion. 

2 Non-extensivity of the average-atom model and 
Tsallis' statistics 

In order to generalize BG statistics, C. Tsallis [12] proposed to replace the 
usual extensive BG entropy 

Sg = -kBY^WilnWi, (3) 

i 

Wi being the probability of the state i of the system, by the entropy 
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1 _ y.w'^ 

ST,q = kB ^^V^, (4) 

q — I 

where q is a positive real number. For instance, turbulence in electron 
plasmas, the flux of solar neutrinos, self-gravitating systems, bare Coulomb 
systems have been successfully described by Tsallis' statistics. BG statistics 
also fails in the interior solar plasma [13] . Given two independent systems in 
the sense of factorizability of the micro-state probabilities, Tsallis' entropy 
of the composite system A + B verifies 



St,M + B) = StAA) + ST,,iB) + (1 - q fTAA)STAB) _ 

Kb 

The quantity |1 — g| quantifies the lack of extensivity of Sx^q (the sys- 
tem is "over-extensive" if g < 1 and "under-extensive" if g > 1). In this 
formalism, the exponential function exp(j;) is replaced by 

exp,(x) = | + (1-^)^1 ^ if [1 + 0; 
I otherwise, 

and the logarithmic function ln(x) by 

ln,(x) = (7) 

The definition of entropjH from Eq. ([5]) warrants that standard BG 
statistics is recovered when g — )• 1. 



3 Determination of the electronic structure 

3.1 Minimization of the non-extensive free energy 

In a non-extensive context, the number of bound electrons is given by 

Nq,b = Y,9k[Uk]', (8) 

k 

where fg^k is a factor characteristic of orbital k, whose degeneracy and 
occupancy are respectively gk and [fg,k]'^- The internal energy reads 

Eg,b = ^9k[fq,k]'' eg^k, (9) 

k 

where €g^k is the energy of orbital k, and the entropy is given by 



^Other g— type distributions exist (see for instance Refs. [141 \T5 
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It is worth mentioning that expression pop comes from the fractal repre- 
sentation of a non-equihbrium electron gas [16j, which states are described 
by a set of orbital occupancy numbers. The total free energy of the plasma 
can be written 



Fq — Fq^b + FqJ + Fq^int + Fq^R 

-^^q(^9^[fq,^r-Z+Z;y (11) 

where Fq^i, = Eq^i, — TSq^ is the free energy of the bound electrons and 
Fq^R the free energy of the radiation field Iq^^- One has 

Fq^R = J hu^'^^^^Yd^dv -T j Sq^^dndu, (12) 

where h is Planck constant and c the speed of light The entropy density 
of the radiation field reads 



q — 1 c-^ 



f\l + ^V-\^V-l\ (13) 



The term Fq^^ represents the free-energy of the delocalized electrons and 
Fq^int the free energy accounting for interactions between electrons. The 
quantity Z* is the average ionization and ^q is the chemical potential, de- 
termined by the preservation of the total number of electrons. Minimization 
of total free energy Fq given by equation (fTT|) with respect to fq^k gives 

[mill] 

= ^ — . (14) 



For q> 1, the distribution ^ gives higher values than Fermi-Dirac dis- 
tribution, and lower values for q <1. Figure [1] represents orbital occupancy 
[fq,kf fo'^ different values of q compared to usual Fermi-Dirac factor f^. The 
bound-electron density reads: 

nq,h{r) = ^gk [fq,kV \'^q,k\^ ■ (1^) 
k 

The minimization of Fq with respect to the delocalized-electron density 
leads to 
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p^dp 



(16) 



where Vq is the new electrostatic potential, 



and 



Xg{r) = yj2m [aq{r) - fiq] 



(17) 



:i8) 



Martinenko and Shivamoggi [18j have shown that non-extensive effects 
on Thomas-Fermi (TF) model can reduce the binding energy, correcting 
in that way a weakness of TF model, due to the divergence r"^/^ of TF 
radial density close to the nucleus, causing a breakdown in the local den- 
sity approximation and a non-physical enhancement of the binding energy 
(boundary effect). However, their expression of the Fermi distribution dif- 
fers from ours: the power is q/{q — 1) instead of l/{q — 1) in our expression. 
Moreover, here we have uncomplete Fermi integrals. Equation ()16p can be 
written 



V2 



[m 



knT 



where 



Fn{a,x,q) 



—dt, 



(19) 



(20) 



'a l + [l + (g-l)(t-x)]9-i 

and the particular integral Fi/2{a,x,q) is convergent only for q < 5/3. 
At low temperatures (x ^> 1), one has 



Fn{x > 1) 



n+1 I 



(l + i)"[l + (g-l)t]^ 
1 + [1 + (^- l)t]^ 



dt. 



which can be written 



FJx > 1) 



X 



n+1 



n + 1 



nx 



n-l 



C2(g) + 



(21) 



(22) 



where Q is defined by 



CM 



ti [l + {q- l)t\ 



■dt. 



(23) 



1 + [1 + (g - l)t]i-^ 

The minimization of the total free energy defined by Eq. (jlip with 
respect to the radiation field Iq^^ gives 
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Figure 3: Ionization versus temperature for a gold plasma at p = 10~^ g/cm^ 
from an LTE calculation and from non-extensive calculations with g=0.98 
and 1.05. 



•■q,v 



1 + 



(g-l) 



hv 



1 

q-l 



(24) 
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the non-extensive radiation field, which does not follow Planck law. Fig- 
ure [2] represents Iq^y for different values of q, compared to usual Planck 
distribution. For q lower (greater) than 1, the maximum is higher (lower) 
than for Planck distribution, and slightly shifted to lower (higher) energy. In 
other words, the energy huQ at which the absorption of photons is maximum 
is larger than 2.822 HbT if g > 1 and lower than 2.822 ksT if g < 1. 
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T 



— LTE, T=20eV 




Ionic charge 



Figure 4: Ionic distributions of an iron plasma at p=10 ^ g/cw? and T=20 
eV from an LTE calculation and from non-extensive calculations with g=0.95 
and 1.05. 
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Figure 5: Ionization for different values of q for an iron plasma at p=\G 
g/cm^ and r=20 eV. 
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3.2 Determination of the self-consistent potential and electro- 
neutrality 



The bound states are obtained solving Schrodinger equation and relativistic 
effects are taken into account in the Pauli approximation [6]. The new 
electrostatic potential Vq can be written: 

Vq{r) = Vg4r) + VqMr), (25) 

where Vq^xc represents the finite-temperature exchange-correlation po- 
tential evaluated following [T2] in the local density approximation (at elec- 
tron density ng{r) = nq^i,{r) + nqj{r)). The coulombic part Vq^c is obtained 
from Poisson equation 

^Vq,c{r) = [nq,b{r)+nqj{r)], (26) 

where n^^fe and Uqj are given respectively by Eqs. (fT5]l and (fT9]) . The 
quantity eo represents the permittivity of vacuum and one has: 

■r->0 

The chemical potential /Ug is determined from the neutrality of the average- 
ion spherical cell: 



7V,,fe = Z-Zl, (28) 

where the average charge Z* is evaluated by integration of the delocalized- 
electron density defined in Eq. (fT9|) . The same process is repeated until 
convergence is reached, i.e. until the potential variation becomes as small 
as required. It is worth mentioning that the use of orbital occupancy fac- 
tors [fq^kf can also be helpful from a numerical point of view [20]. Indeed, 
average-atom models can suffer from convergence problems at low tempera- 
ture, due to the fact that the Fermi-Dirac distribution tends to a Heaviside 
function. Such difficulties might be avoided performing the self-consistent 
calculation with [fq^kf distribution, and taking the limit q ^ 1. 



4 Orbital energies, populations and average ion- 
ization 

Tables [1] and [2] display respectively the energy and population of several 
orbitals for an iron plasma at /9=10^^ g/cm^ and T=20 eV from an LTE 
calculation and from non-extensive calculations with g=0.95 and 1.05. When 
q <1 (respectively g > 1), energies and populations are higher (respectively 
lower) than in the LTE case. Figure [3] represents the average ionization 
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Orbital 


Energy 


Energy 


Energy 




at LTE 


for g=0.95 


g=1.05 


3s 


-238.90 


-205.71 


-289.51 


3Pl/2 


-204.89 


-171.97 


-255.77 


31^3/2 


-146.66 


-114.01 


-198.36 


4^5/2 


-60.95 


-43.54 


-88.98 


6/5/2 


-16.60 


-10.07 


-27.75 


8/I11/2 


-5.99 


-2.99 


-11.16 



Table 1: Energies (eV) of a few selected orbitals for different approaches in 
the case of an iron plasma at p=10^^ g/cm^ and T=20 eV: LTE case and 
non-extensive cases with q=0.95 and g=1.05. 

(1995) 



Orbital 


Population 


Population 


Population 




at LTE 


for q=0.95 


for q=1.05 


3s 


1.92 


1.93 


1.84 


3pi/2 


1.62 


1.69 


1.38 


3d3/2 


0.77 


0.92 


0.45 


4d5/2 


0.02 


0.05 


0.003 


6f5/2 


0.002 


0.0099 


0.0001 


8hii/2 


0.002 


0.014 


0.0001 



Table 2: Population (number of electrons) of a few selected orbitals for 
different approaches in the case of an iron plasma at p=10~^ g/cm^ and 
r=20 eV: LTE case and non-extensive cases with q=0.95 and g=1.05. 

versus temperature for a gold plasma at p=10~^ g/cm^ in LTE and in the 
non-extensive formalism for §=0.98 and 1.02. We can see that the non- 
extensive ionization is lower than the extensive one for q < 1 and larger 
for q > 1. Figure S] represents the ionic distributions in an iron plasma 
at p=10-^ g/cm^ and r=20 eV in the LTE case and calculated with the 
gi— version of the code for q=1.05 and g=0.95. The non-extensive distribu- 
tions are shifted to lower charge states {if q < 1) and higher (if q > 1) and 
their asymetry can be different from the LTE one: the ionic distribution 
is wider for q < 1. The difference increases with the value of q. Figure [5] 
shows non-extensive ionization Z* versus non-extensivity parameter q. One 
can notice two steps: the first one corresponds to Z*=16, i.e. 10 remain- 
ing bound electrons, corresponding to full K and L shells. The second step 
corresponds to the asymptote of full ionization Z* — )• Z=26. 
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5 Conclusion 



In order to investigate the impact of non-extensivity on atomic physics in 
hot plasmas, a self-consistent atomic-structure model was presented, result- 
ing from the minimization of a non-extensive free-energy and including the 
radiation field. All electrons (bound and delocalized) are described in the 
framework of the non-extensive formalism. Delocalized electrons are taken 
into account using non-extensive semi-classical Thomas-Fermi approxima- 
tion, which leads to a reduction of the binding energy, correcting in that way 
a weakness of standard Thomas- Fermi model. Standard thermodynamics is 
recovered when g — )• 1 and the formalism presented here can be extended to 
other approaches of hot-plasma atomic structure. For instance, another pos- 
sibility to account for the electron states consists in considering not only the 
mean configuration of the plasma, but the real configurations (with integer 
electron populations for the orbitals). Because the number of configurations 
is huge, they can be gathered into superconfigurations [211 1221 1231 [Ml [25] 
and the self-consistent calculation of electronic structure presented in this 
paper can be performed for each superconfiguration. In view of the growing 
interest for Tsallis' non-extensive statistics, we hope that the present study 
will give rise to new ideas for the simulation of plasma atomic physics. 
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